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FAST  IONS  FROM  LASER  PLASMAS : ANALYTIC  SOLUTION  AND  SCALING  LAWS 

Fast  ions  are  a characteristic  phenomenon  observed  in  Irradiation 
of  targets  by  high-power  lasers.  Electrostatic  acceleration  by  a steep, 
localized  pressure  gradient  with  a high  electron  temperature  is  a 
suggested  mechanism  for  producing  the  fast  ions^”^.  Recent  work  at 
NRL^  demonstrates  that  a multi-species,  amblpolar  expansion  model  can 
account  for  most  of  the  high-energy  ion  distribution  measured  far  from 
a Nd-laser  target.  However,  the  ions  of  highest  energy  can  be  treated 
more  effectively  with  the  analytic  solution  presented  in  this  note. 

The  only  previous  analytic  solution  describing  plasma  expansion  by 
this  process  of  electrostatic  acceleration  is  a self-similar  solution'^. 
Neither  the  self-similar  solution  nor  the  closely  related  numerical 
solution^  involves  a characteristic  length,  and  neither  yields  scaling 
laws  for  the  maximum  ion  energies,  since  both  feature  ion  velocities 


that  increase  indefinitely. 

In  this  note  we  present  a new  analytic  solution  that,  unlike  the 
self-similar  solution,  involves  a characteristic  scale  height,  and 
yields  an  expression  for  the  maximum  energy  of  ion  species  in  terms  of  . — 
their  charge,  mass,  and  scale  height,  and  the  electron  temperature.  '"'Court 

o'- 

We  also  suggest  a means  by  which  the  expansion  described  by  the  new 
solution  might  evolve  into  a self-similar  expansion,  - — 

Note:  Manuscript  submitted  August  10,  1977.  ! ORUlB'MK'AVAWBilin  CO^ 


1 


Oist.  A 

V.!i.  J.;J/'or 

A 

□ □ 


Our  assumpClons  and  equations  describing  the  expanding  plasma  are 
essentially  those  of  Refs.  5 and  '•*.  We  begin  by  describing  the  plasma 
by  a two-fluid  model  involving  hot  electrons  and  relatively  cold  ions. 
Later  we  generalize  the  analytic  solution  to  more  than  one  ion  species. 
Each  fluid  satisfies  a continuity  equation.  The  quasi-neutrality 
condition  replaces  the  Poisson  equation,  because  the  characteristic 
scale  height  is  taken  to  be  much  greater  than  the  electron  Debye  length. 
The  electrons  are  assumed  to  have  a Boltzmann  distribution  (electron 
pressure  balanced  by  amblpolar  potential)  with  constant  temperature. 

The  validity  of  assuming  isothermal  electrons  depends  upon  high  thermal 
conductivity  from  a heat  reservoir  supplied  by  the  laser.  The  electron 
temperature  is  assumed  high  enough  (~  tens  of  keV)  to  neglect  electron- 
ion  collisions.  The  ion  temperature  is  assumed  to  be  high  enough  to 
neglect  viscosity,  but  much  less  than  the  electron  temperature,  so  that 
the  ion  pressure  can  be  ignored.  Consequently,  the  ion  momentum 
equation  has  only  the  ambipolar  electric  field,  or  equivalently  the 
electron  pressure  gradient,  as  its  source  term. 

The  system  of  equations  in  plane  geometry  reduces  to 

^(In  n)/at  + V a(-tn  n)/ax  + c 5(v/c)/Sx  = 0, 

5(v/c)/at  + V a(v/c)/ax  + c a(-tn  n)/ax  = o,  (i) 

in  which  n is  the  ion  density  and  v the  velocity  of  the  ion  fluid. 

The  constant  sound  speed  is  c s (ZT^/M)^,  in  which  Z and  M are  the 
charge  state  and  mass  of  the  ions,  and  T^  is  the  electron  temperature 
in  energy  units.  These  equations  can  be  posed  as  a Cauchy  initial-value 
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problem,  and  solved  generally  by  the  method  of  characteristics,  but 
analytic  solutions  are  useful. 

Our  initial  conditions  are  n = n^  exp(-x/^.),  and  v = C,  in  which 
n^  and  are  a constant  density  and  scale  height.  These  initial 
conditions  are  consistent  with  a quick  deposition  of  some  laser  energy 
at  the  target  surface  (during  the  rise  of  a laser  pulse  for  example). 
With  these  initial  conditions,  our  analytic  solution  to  the  eqs.  (1)  is 

n = n^  exp(-x/.t  + v = (2) 

In  a reference  frame  having  acceleration  c^/(,  , the  solution  be- 
comes n = n^  exp(-x/^,)  and  v = 0.  Thus  the  analytic  solution  describes 
an  exponential  density  profile  maintaining  its  shape  as  it  accelerates 
at  a rate  c^lK.  In  the  lab  frame  the  energy  per  ion  is 

E = 4Mv^  a i(Z^/M)(T  /'t)^t^  (5) 

e 

Here  is  evident  the  2^/M  scaling  law  found  elsewhere^  theoretically 
and  experimentally.  Moreover  eq.  (3)  demonstrates  that  fast  ions  are 
produced  by  high  electron  temperatures  and  steep  density  gradients.  An 
estimate  for  the  acceleration  time  scale  (and  therefore  the  peak  energy) 
will  be  provided  below. 

If  the  initial  conditions  in  some  region  of  positive  x are  Instead 
n = n^  exp(-x/^)  and  v = cx/t,  then  the  solution  to  eqs.  (1)  is  the 
self-similar  solution*^  n = n^  expr-(x  + ct)/(.?  + ct)]  and 
V = c(x  + ct)/(-t  + ct).  However  if  the  density  gradient  extends  over 
several  scale  heights,  then  these  initial  conditions  imply  initial 
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ion  velocities  of  several  times  sound  speed,  or  ions  that  are  already 
fast.  Thus,  the  initial  conditions  for  the  self-similar  solution 
require  some  strong  acceleration  mechanism  to  have  been  in  operation 
prior  to  t = C,  while  the  initial  conditions  leading  to  the  solution 
(2)  do  not.  We  now  provide  a hybrid  analytic-numerical  example  to 
illustrate  how  the  solution  (2)  can  provide  the  acceleration  mechanism 
necessary  to  launch  the  self-similar  plasma  expansion. 

Let  two  families  of  curves,  and  , be  defined  in  the  x,t 
plane  by  dx/dt  = v + c and  dx/dt  = v - c respectively.  are  the 
forward -running  characteristics,  and  Cg  the  backward-running  charac- 
teristics. The  Riemann  invariants  are  v + c-tn(n/n^) , which  has  a 
constant  value  along  C^,  and  -v  + c/,n(n/n^),  which  has  a constant  value 
along  Cg.  If  the  values  of  n and  v are  specified  along  a segment  of 
the  line  t = C in  the  x,t  plane,  then  with  the  aid  of  the  Riemann 
invariants,  (1)  can  be  solved  as  a Cauchy  initial-value  problem®. 

For  example  consider  a density  profile,  initially  at  rest,  given 
by  n = n^  exp(-x/A)  for  x s 0,  and  n = n^  for  x « 0,  as  illustrated  by 
the  dashed  line  in  Fig.  la.  The  characteristics  are  shown  in  Fig.  lb 
in  the  two  regions  where  analytic  solutions  are  known.  In  region  I 
the  characteristics  are  straight  lines,  indicating  that  the  plasma 
remains  at  rest  at  constant  density  there.  In  region  III  the  equations 
of  the  and  Cg  characteristics  are  x = c^t^/2l^  + ct  + and 
X = c^t^/2/,  -ct  + Xg  respectively,  in  which  x^  and  x^  are  constants. 

In  this  region  the  fluid  accelerates  to  the  right  at  a rate  ! K , and 
the  solution  (2)  pertains. 


At  t = 0,  the  signal  communicating  the  abrupt  change  in  density 
at  X = C propagates  to  the  left  at  sound  speed  c,  and  to  the  right  at 
sound  plus  flow  speed,  c + z^tlK  . Within  the  expanding  region  II, 
bounded  by  these  two  signal  fronts,  the  solution  can  be  seen  in  Fig. 
la  to  be  nearly  self-similar  (linear  velocity  and  logarithmic  density 
profiles) , 

The  solution  in  Fig.  la  was  calculated  by  a 1-D  hydrodynamics 
code^  from  the  initial  profiles  shown  by  the  dashed  curves.  In  the 
region  to  the  right  of  the  signal  emanating  from  x = 0,  the  velocity 
plateau  calculated  by  the  code  agrees  exactly  with  the  analytic  value 
given  by  (2),  namely  v = This  example  therefore  suggests  that 

the  solution  (2)  describes  a mechanism  whereby  ions  can  be  accelerated 
from  rest  to  the  velocities  necessary  to  initiate  the  self-similar 
expansion  of  a plasma. 

In  a reference  frame  having  acceleration  c^/^t  , the  fluid  in  region 
III  would  appear  to  be  at  rest,  and  the  and  would  be  straight 
lines  with  dx/dt  = ±c.  Therefore  any  localized,  exponential  density 
profile  of  width  L is  eroded  at  each  end,  and  extinguished  after  a 
time  L/2c.  For  example  in  Fig.  lb  it  can  be  seen  that  an  exponential 
profile  of  width  L = is  extinguished  after  a time  2t/c.  Abrupt 

changes  in  the  profiles  in  Fig.  la  appear  at  the  locations  of  the 
signals  that  originated  at  x = C. 

If  more  than  one  ion  species  is  present,  and  each  is  treated 
within  the  context  of  the  fluid  model  described  earlier,  then  their 
momentum  equations  decouple  in  any  region  in  which  each  species  starts 
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from  rest  with  an  exponential  density  profile  of  the  same  scale  \ 

height.  In  such  a region  each  ion  species  behaves  as  though  the 

others  were  not  present,  and  the  solution  (2)  is  valid  for  each  species.  ! 

If  such  a region  has  width  L,  then  it  is  extinguished  after  a ^ j 

time  tc.  , = L/2c  , in  which  c is  the  sound  speed  of  the  multi-  , 1 

final  s’  s f I 

species  fluid.  Substituting  in  (3),  we  obtain  a lower  bound 

for  the  peak  ion  energy  of  any  ion  species,  assuming  that  each  species 
is  accelerated  from  coexisting  exponential  density  profiles  with  scale 
height  I , 


Thus  if  the  exponential  density  profile  is  a ‘few  scale  heights  wide, 
the  directed  energies  of  the  ions  can  exceed  the  electron  thermal 
energy  before  the  self-similar  expansion  becomes  the  dominant  accelera- 
tion process.  However  if  the  ion  velocity  becomes  comparable  with  the 
electron  thermal  velocity,  then  the  assumption  of  a Boltzmann  distri- 
bution breaks  down,  and  this  model  is  no  longer  applicable. 

In  conclusion  we  have  found  an  analytic  solution  for  the  ambipolar 
expansion  of  a cold,  collisionless  ion  fluid  from  a steep  pressure 
gradient  with  a hot  isothermal  electron  background.  The  solution,  in 
contrast  to  the  self-similar  solution,  describes  a means  for  accelerating 
ions  from  rest  while  maintaining  their  exponential-density  scale  height. 
However,  a steep,  localized,  exponential  gradient  is  eroded  at  both 
ends  at  sound  speed,  and  finally  extinguished,  probably  leaving  behind 
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a plasma  accelerating  in  an  approximately  self-similar  expansion,  A 
scaling  law  for  the  peak  ion  energies  in  a multi-species  ion  fluid 
was  derived  that  shows  a strong  dependence  on  electron  temperature  and 
density  scale  height,  as  well  as  the  Z^/M  scaling. 

The  authors  wish  to  express  their  gratitude  to  S.  E,  Bodner  and 
B,  H,  Ripin  for  helpful  discussions. 
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Fig.  1 — (a)  Velocity  (in  units  of  sound  speed  c)  and  logarithmic  densit>’ 
profiles  initially  (dashed)  and  at  time  2ii/c  (solid),  C is  scale  height,  and 
(b)  characteristic  curves  for  same  initial  conditions  as  (a).  Analytic  solu- 
tions of  density  and  velocity  in  regions  I and  III  are  given  in  text. 
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